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ABSTRACT 

Observations and numerical magnctohydrodynamic (MHD) simulations indicate the 
existence of outflows and ordered large-scale magnetic fields in the inner region of 
hot accretion flows. In this paper we present the self-similar solutions for advection- 
dominated accretion flows (ADAFs) with outflows and ordered magnetic fields. Stimu- 
lated by numerical simulations, we assume that the magnetic field has a strong toroidal 
component and a vertical component in addition to a stochastic component. We obtain 
the self-similar solutions to the equations describing the magnetized ADAFs, taking 
into account the dynamical effects of the outflow. We compare the results with the 
canonical ADAFs and find that the dynamical properties of ADAFs such as radial 
velocity, angular velocity and temperature can be significantly changed in the pres- 
ence of ordered magnetic fields and outflows. The stronger the magnetic field is, the 
lower the temperature of the accretion flow will be, and the faster the flow rotates. 
The relevance to observations is briefly discussed. 

Key words: accretion, accretion discs - magnetohydrodynamics: MHD ~ ISM: jets 
and outflow - black hole physics 



1 INTRODUCTION 

Advection-dominated accretion flows (hereafter ADAFs) 
have been studied extensively (e.g, Narayan & Yi 1994, 
1995; Abramowicz et al. 1995; see Narayan, Mahadevan & 
Quataert 1998 and Kato, Fukue & Mineshige 1998 for re- 
views). It is now rather well established that this accretion 
mode exists in the quiescent and hard states of black hole 
X-ray binaries and low-luminosity active galactic nuclei (see 
Narayan 2005, Yuan 2007, Ho 2008, and Narayan & Mc- 
Clintock 2008 for recent reviews). ADAFs only exist below 
a critical mass accretion rate. Below this accretion rate, the 
rate of radiative cooling is so weak that the viscous heating 
is balanced by the advection. With the increase of accretion 
rate, radiation becomes more and more important until it 
becomes equal to the viscous heating at this critical rate. In 
this case, the energy advection is equal to zero (Narayan, 
Mahadevan & Quataert 1998). Above this critical accretion 
rate, up to another limit close to the Eddington accretion 
rate, another hot accretion solution — luminous hot accretion 
flows (LHAFs) — was found which is a natural extension of 
ADAFs (Yuan 2001, 2003). In this solution, the flow is able 
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to remain hot because the radiative cooling rate, although it 
is strong, is still lower than the sum of compression work and 
viscous heating. Note that the cool thin solution — standard 
thin disk — always exists in the whole range of accretion rates 
of ADAFs and LHAFs. In the present work we only focus on 
the hot accretion flows — ADAFs, but our discussion should 
hold for LHAFs. In the early version of ADAFs, the accre- 
tion rate is assumed to be a constant, i.e., there is no outflow 
throughout the region. The magnetic field is included, but 
only its stochastic component, while the large-scale ordered 
field is not considered. 

Great developments have been achieved since the dis- 
covery of ADAFs. One is the realization of the existence 
of outflows in ADAFs, i.e, only some fraction of the accre- 
tion material available at the outer boundary is actually ac- 
creted into the central black hole (Narayan & Yi 1994, 1995; 
Blandford & Begelman 1999; Stone, Begelman, Pringle 1999; 
Igumenshchev & Abramowicz 1999; Stone & Pringle 2001). 
The physical reason for the origin of the outflow is believed 
to be that the Bernoulli parameter of the flow is positive 
(Narayan & Yi 1994; Blandford & Begelman 1999). Another 
possible mechanism of the outflow origin is associated with 
large-scale ordered magnetic fields, e.g., through the magne- 
tocentrifugal force (Blandford & Payne 1982; Henriksen & 
Valls-Gabaud 1994; Fiege & Henriksen 1996). Observations 
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in the center of our Galaxy supply strong evidence for the ex- 
istence of outflows. From Chandra observations combined 
with Bondi accretion theory, we can predict the accretion 
rate at the Bondi radius. Polarization observations at ra- 
dio wavebands, however, indicate that the accretion rate at 
the innermost region must be significantly smaller than the 
Bondi value (Yuan, Quataert & Narayan 2003). Therefore, 
a large amount of material must be lost into outflows. 

Another interesting result of numerical magnetohydro- 
dynamic (MHD) simulations of the hot accretion flow is that 
a large-scale ordered magnetic field exists in the inner re- 
gions of ADAFs. Independent of the initial configuration 
of the magnetic field (toroidal or poloidal) in the main 
body of the accretion flow the field is primarily toroidal, with 
weak radial and vertical components. This large-scale struc- 
ture is imposed on the stochastic component of the magnetic 
field on small scales (Machida, Hayashi & Matsumoto 2000; 
Hirose et al. 2004). 

Both outfiows and large-scale magnetic fields can af- 
fect the dynamics of ADAFs significantly. For example, both 
of them can effectively transfer angular momentum. These 
are alternative mechanisms in addition to the turbulence 
mechanism associated with the magnetorotational instabil- 
ity (MRI) proposed by Balbus & Hawley (1991; 1998). Stone 
& Norman (1994; see also Mouschovias & Paleologou 1980) 
investigate the angular momentum transfer by magnetic 
braking effects associated with a large-scale magnetic field. 
If the specific internal energy of the outflow is different from 
that of the inflow where the outflow originates, the outflow 
acts as an extra cooling or heating term in the accretion flow, 
as discussed phenomenologically by Blandford & Begelman 
(1999). Xie & Yuan (2008) parameterize the outflow proper- 
ties and systematically investigate the effects of the outflow 
on the dynamics of the inflow, in absence of the large-scale 
magnetic field. 

It is thus necessary to investigate the dynamics of 
ADAFs with coexistent outflows and large-scale magnetic 
fields. Several works have been done recently under the 
self-similar approximation (Akizuki & Fukue 2006; Abbassi, 
Ghanbari & Najjar 2008; Zhang & Dai 2008) or global solu- 
tion (Oda et al. 2007). All these works consider the dynam- 
ical effects of the outflow by adopting the form M oc r"^^^'^ 
(e.g., Eq. ([2TJ; Blandford & Begelman 1999) to describe the 
accretion rate while all other effects such as the probable an- 
gular momentum transfer by outflows are neglected. In Ak- 
izuki & Fukue (2006), Oda et al. (2007), and Abbassi, Ghan- 
bari & Najjar (2008), only the toroidal component of the 
large-scale magnetic fleld is considered; thus, the large-scale 
magnetic field in their model only supplies an additional 
force in the radial direction, while it is unable to transfer 
angular momentum. In Zhang & Dai (2008), although all 
the three components of the large-scale magnetic fleld are 
included explicitly, their solutions unfortunately violate the 
magnetic divergence-free condition when s 7^ 0. 

In this paper, we investigate the self-similar solutions 
of ADAFs with coexistent outflows and large-scale magnetic 
fields. We assume the large-scale magnetic field has both z 
and 4> components, hence, the large-scale magnetic fleld does 
not only affect the radial force balance of the accretion flow, 
but also helps to remove angular momentum. Following Xie 
& Yuan (2008), we take into account the dynamical effects 
of outflows by considering the differences of specific angular 



momentum and specific internal energy between inflow and 
outflow. The paper is organized as follows. In §2 we present 
the basic MHD equations, which include ordered magnetic 
fields and outflows. We address in §3 further assumptions, 
and achieve a set of self-similar equations. Numerical results 
are presented in §4. We summarize the paper in §5. 



2 BASIC EQUATIONS 

In cylindrical coordinates (r, (f>, z), we investigate the steady- 
state, axisymmetric [d/dt — d/d4> — 0) magnetized 
advection-dominated accretion flows with outflows. The gen- 
eral MHD equations in Gaussian units read (e.g., Lovelace, 
Romanova, & Newman 1994), 
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The induction equation of the magnetic fleld reads, 
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The final equation is the energy equation. 
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Here is the gravitational potential; Fvis is the viscous force 
per unit volume; B is the magnetic field; is the magnetic 
diffusivity; J is the current density; e is the specific internal 
energy of the accretion flow; q^is and Qrad are the viscous 
heating and radiative cooling rates; p, v and p have their 
usual meanings. For simplicity, we use the Newtonian po- 
tential and neglect the self- gravity of the accreting gas. The 
advection factor, / (0 ^ / ^ 1), describes the fraction of the 
viscous dissipation energy which is stored in the accretion 
flow and advected into the central black hole rather than 
radiated away. 

From the numerical MHD simulations (e.g., Machida, 
Hayashi & Matsumoto 2000; Hirose et al. 2004), it is rea- 
sonable to decompose the magnetic fleld into one large-scale 
component and one fluctuating or turbulent component. The 
general picture that emerged from the simulations is that the 
main body of the accretion flow is governed by a toroidal 
magnetic field, especially in its inner region, while regions 
near the poles are primarily governed by a poloidal mag- 
netic field (see Fig. 6 in Hirose et al. 2004), which is mainly 
in the vertical direction. Thus, we only consider the and 
z components [3^,3^) of the large-scale magnetic field and 
neglect the radial component (-Br). We describe the effects 
of the fluctuating component of the magnetic fleld in trans- 
ferring the angular momentum and in dissipating the energy 
through the usual a description. Speciflcally, the viscous 
heating rate q^is in Eq. (|6]) is associated with the turbulent 
component instead of the large-scale ordered component of 
the magnetic field. This is because numerical simulations 
show that the dissipation mainly comes from the thermal- 
ization of the magnetic energy via magnetic reconnection 
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(Hirose et al. 2006), which is obviously associated with the 
turbulent component. 

Following Lovelace, Romanova & Newman (1994), we 
assume is an even function of z. Furthermore, we neglect 
the vertical gradient of in the disc. We take to be an 
odd function of z. We therefore have, 

f^^O, (7) 
oz 



B^ — Bo — , B^\z=H ~ —B^\z=^H , 

£1 

Bo ~ B^\z^H, 



(8) 
(9) 



where H is the half-height of the accretion disc defined as 
H = Csr/vK{l + /3i)^''^, where Cs is the isothermal sound 
speed, vk is the Keplerian velocity and /3i is defined below in 
Eq. (22). It will reduce to its traditional form of H = Car/vK 
in absence of the toroidal component (B^) of the ordered 
magnetic field (Eq. (|13|) ). Under the above assumptions, the 
divergence- free condition of the magnetic field (Eq. (j4|) is 
automatically satisfied. 

The specific angular momentum and specific internal 
energy of the outflow are in general different from those of 
the inflow where the outflow originates. A two-dimensional 
calculation is ideal but complicated and out of the scope of 
our current work. Instead, we follow the method illustrated 
in Xie & Yuan (2008) and use several parameters (^i,^2,C3i 
see below for definitions) to characterize the properties of 
the outflow. In this way, we integrate the above MHD equa- 
tions in the vertical direction and get the following 1.5- 
dimensional equations to describe the inflow dynamics (Xie 
& Yuan 2008). 

The equations of the conservation of mass and momen- 
tum are, 



dMjr) 
dr 



Here Vr,w, 1^0, m and Vz,in denote the r, (j) and z components 
of the velocity of the outflow when it is just launched from 
the accretion flow (Xie & Yuan 2008), and e-w and e corre- 
spond to the speciflc internal energy of the outflow and that 
of the inflow, respectively. Following Xie & Yuan (2008), we 
assume Vr,-w = ^ifr,i'0,ii, = ^2V^ and = ^35. Parameter 
7 is the ratio of the speciflc heats, which goes from the rel- 
ativistic value 4/3 to the typical monatomic gas value 5/3. 
In our calculations we set 7 = 4/3. We flnd our results are 
not sensitive to the value of 7. 

The second terms on the left-hand side of Eqs. (Illf) & 
(|12|) are the momenta (radial and angular) taken away from 
or deposited into the inflow by the outflow (Xie & Yuan 
2008), and the third term on the left-hand side of Eq. (fT4t 
is the internal energy taken away from or deposited into the 
inflow by the outflow. If all these ^s are equal to unity, our 
model will reduce to those neglecting the discrepancies of 
speciflc momenta and internal energy between the inflow and 
the outflow (e.g., Akizuki & Fukue 2006; Abbassi, Ghanbari 
& Najjar 2008; Zhang & Dai 2008). 

Now we discuss the induction equation (Eq. Q). This 
equation describes the growth or escape rate of the magnetic 
fleld by dynamo and diffusion. We deflne the advection rates 
of the two components of the magnetic flux as, 
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For a steady-state accretion system considered here, if the 
dynamo and the diffusion effects are neglected, the above 
two quantities will be constant. However, in a realistic ac- 
cretion disc, the above two quantities can vary with radius 
due to the presence of the dynamo and the diffusion effects. 
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The static equilibrium in the vertical direction is achieved 
between the ^-component of the gravitational force and the 
(gas and magnetic) pressure gradient force. 
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The tension term associated with B^ is BrdB^/dr. It is ne- 
glected since we have assumed Br = 0. The energy equation 
is 

Vr del "rCs dp 1 dM (r) 
7 — 1 dr p dr 27rrE dr ^ 
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3 SELF-SIMILAR SOLUTIONS 

We seek self-similar solutions in the following forms (e.g.. 
Narayan & Yi 1995), 



W0 



-ci avK, 



C2 Vk, 



Cs = C3 Vk- 



(17) 
(18) 
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We assume the surface density (E = 2pH) to have the form, 

E = Eo r^ (20) 

where Eo and s are two constants, so the accretion rate takes 
the form, 



M 



(21) 



Note that s = —1/2 will lead to solutions without outflows. 

We assume the magnetic pressure on the surface of the 
accretion flow to be proportional to the gas pressure in the 
following way. 



{B^\z=Hf/8n 
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Pi and (52 are two parameters. The induction equations will 
be satisfied if tlie dynamo and tfie diffusion of tlie magnetic 
field satisfy the self-similar conditions required by the ad- 
vection rate (Eqs. (|15|l & (|16|) ). Typical values of /3i and 
j32 lie in the range 0.01 — 1 in hot accretion flows (e.g., 
De Vilhers, Hawley & Krohk 2003; Beckwith et al. 2008), 
but here we also consider the magnetically-dominated case ( 
/3i, /?2 > 1). On one hand, the MHD numerical simulation by 
Machida, Nakamura & Matsumoto (2006) shows that when 
thermal instability happens in an ADAF, the thermal pres- 
sure rapidly decreases while the magnetic pressure increases 
due to the conservation of magnetic flux. This will result in 
large /3i and P2 and forms a magnetically-dominated accre- 
tion flow. In addition, an initially weak large-scale fleld at 
the outer region can be amplifled due to the compression 
of the flow as it is accreted inward. This can also result in 
large fii and /32. Of course, in this case, the MRI will be 
suppressed, but we still keep a in equations for simplicity 
and this treatment will not affect our results significantly 
for our purpose. 

The momentum conservation and the vertical equilib- 
rium equations (Eqs. (|ll| - [T3t l now reduce to, 
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Obviously for given values of a,f,s,f5\ and (32, Eqs. pi)) . 
1(25)1 and (f27)l form a closed set of equations of ci , C2 and C3 , 
which will determine the dynamics of the accretion flow. 



4 RESULTS 

The first and the second terms on the right-hand side of 
Eq. (|25p are angular momentum transfer due to turbulent- 
viscosity (MRI) and the large-scale magnetic field, respec- 
tively. We define a new parameter C4 to account for the ratio 
of angular momentum transport due to these two mecha- 



C4 



3a(H-s)Vl + /3iC2^ 



(28) 



The large-scale magnetic fleld will be the dominant mech- 
anism in angular momentum transport when C4 is below 
unity. 

We flrst show in Fig.[T]the effects of the large-scale mag- 
netic fleld on the accretion flow. We set a = 0.1, / = 1.0 and 
s = —0.5 (no outflows). The solid, dotted and long-dashed 
lines correspond to /32 = 0.01, 1.0 and 10, respectively. We 



can see from the flgure that as {/3i) increases, the radial 
(ci) and angular (02) velocities increase, while the temper- 
ature (cs) and the turbulent-viscosity contribution to an- 
gular momentum transport (04) decrease. The decrease of 
C4 is obvious and easy to understand, since the large-scale 
magnetic fleld should be more important in transporting an- 
gular momentum as its strength increases. We flnd that the 
magnetic pressure gradient force (the fourth and flfth terms 
on the right-hand side of Eq. (|24|l ') is a centrifugal force, 
while the magnetic stress force (the sixth term on the right- 
hand side of Eq. (|24|) ') is a centripetal force. For the wide 
range of /3i and P2 considered here, the gradient of the total 
pressure (gas plus magnetic pressure) serves as a centrifu- 
gal force, and it exceeds the magnetic centripetal force. All 
these terms are proportional to the sound speed (or, equiva- 
lently, the temperature C3). The solutions to Eqs. (|24|) - (|27p . 
as shown in Fig. [T] indicate that an increase in {j3i ) re- 
sults in a decrease in temperature (03). This in turn leads 
to a decrease of the effective centrifugal force, and subse- 
quently an increase of the infalling velocity Vr, as shown in 
the upper-left panel of Fig. \T\ Combined with the energy 
equation (Eq. (|27|) ). the angular velocity will increase corre- 
spondingly. 

The change of the z component of the ordered mag- 
netic fleld has similar effects. It is evident from Figure [1] 
that an increase in the z component of the magnetic fleld 
Bz {P2) will lead to a decrease in temperature (C3). Due to 
the decrease of temperature, the effective centrifugal force 
decreases, so the infalling velocity increases. From the en- 
ergy equation (Eq. (|27p ). the angular velocity will increase 
correspondingly. We thus come to a general conclusion that, 
the stronger the magnetic fleld is, the faster the accretion 
flow rotates and falls into the central object, and the lower 
the temperature will be (see Figs. ((2} & (|3]) for situations 
with outflows). 

The decrease of temperature of ADAFs after taking 
into account the large-scale magnetic fleld may have an im- 
portant observational implication. The X-ray spectra of the 
hard state of black hole X-ray binaries and some AGNs (such 
as type 1 Seyfert galaxies) are well described by a power- 
law form with a high-energy cutoff. This kind of spectrum 
is widely accepted as originating from the thermal Comp- 
tonization process in a hot plasma which is likely to be the 
region of a hot accretion flow where most of the radiation 
originates (see review by Zdziarski & Gierlinski 2004). The 
value of the cutoff energy is determined by the temperature 
of the accretion flow, while the slope of the power-law dis- 
tribution is determined by the product of temperature and 
optical depth. From fltting the observed X-ray spectrum, we 
can well constrain these two parameters, namely the opti- 
cal depth and temperature. On the other hand, we can also 
calculate their values from the hot accretion flow model. 
Yuan & Zdziarski (2004) flnd that while the theoretical pre- 
diction is roughly consistent with observational results, for 
some sources the temperature obtained from observation, 
~ 40 — 60 keV, is signiflcantly lower than that the lowest 
value predicted by the hot accretion flow models, > 80 keV 
(see Fig. 1 in Yuan & Zdziarski 2004). Our current work 
indicates that this discrepancy may be solved by including 
the large-scale magnetic fleld (Fig.[T| and the outflow (Figs. 
[2]y3|. Unfortunately, it is difficult to give a quantitative es- 
timation of how strong the magnetic field is required to be 
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Figure 1. The dynamics of tlie accretion flow in absence of outflows. Parameters are set as s = — 0.5,o = 0.1 and / = 1. The solid, 
dotted and long-dashed lines correspond to /32 = 0-01, 1 and 10, respectively. The bottom-right panel shows the ratio of the angular 
momentum transport due to turbulence (MRI), to that due to the large-scale magnetic field (see Eq. II28I I for definition). 



to solve the discrepancy. This is because our current work is 
based on the self-similar assumption which does not hold in 
the innermost transonic region of the accretion flow where 
most of the radiation originates. 

Before examining in detail the effects of the outflow on 
the dynamics of the inflow, we first check the effects of the 
discrepancy of the radial velocity between the outflow and 
the inflow (indicated by ^i). Consistent with Xie & Yuan 
(2008), we find that for relatively weak outflows (e.g., s < 0) 
and reasonable values of ^i, the effects are of minor impor- 
tance. Thus for the discussions below, we set = 1. 

We present in Fig. [5] the effects of angular momentum 
discrepancy between the inflow and the outflow. We set q = 
0.1,/ = = 6 = 1 and fi-z = 0.01. Parameter s = 0.5 
means Af(r) oc r, which corresponds to a relatively strong 
outflow (Blandford & Begelman 1999). The dotted, solid 
and long-dashed lines correspond to ^2 = 0.8, 1 and 1.2, 



respectively. Parameter ^2 > 1 means that the outfiow helps 
to remove angular momentum from the inflow. 

The accretion flow has distinctive properties between 
Pi < 1 and /3i > 1 . When /3i is below unity, the effects of the 
large-scale magnetic field can be neglected. As ^2 increases 
from 0.8 to 1.2, we can see from Fig. [2] that the radial and 
angular velocities increase while the temperature decreases. 
We can understand these results in the following way. From 
Eq. (|27|l . we see that ci and C2 will increase or decrease to- 
gether (note the magnetic terms are negligible). If both of 
them would decrease, from Eq. (24) we know that C3 must 
increase; while this is in conflict with what Eq. (25) implies 
if ^2 increases. Thus both ci and C2, i.e., the radial and an- 
gular velocities, must increase when ^2 increases, and the 
temperature decreases. Physically, the increase of the radial 
velocity is explained by the decrease in the gas pressure gra- 
dient force with the decreasing temperature. We note that 
these self-similar results are consistent with the global solu- 



6 D.F. Bu, F. Yuan, and E.G. Xie 



-0.2r 



. 4r 



-0.4- i- 0.3 




-2 -1.5 -1 -0.5 0.5 1 -2 -1.5 -1 -0.5 0.5 1 



Figure 2. The dynamics of tlie accretion flow when the outflow has different angular velocity from that of the inflow. Parameters are 
/32 = 0.01, a = 0.1, / = 1, s = 0.5 and ^1=^3 = 1.0. The accretion rate M is then proportional to r. The dotted, solid and long-dashed 
lines correspond to §2 = 0.8,1 and 1.2, respectively. The bottom-right panel shows the ratio (C4; Eq. ||28| |) of the angular momentum 
transport due to MRI, to that due to the ordered magnetic field. 



tions of Xie & Yuan (2008). We can understand the results 
for /3i > 1 in a similar way, except that in this case the mag- 
netic force (magnetic pressure and stress forces) exceeds the 
gas pressure gradient force. The critical value of /3i corre- 
sponds to the balance between the forces of gas pressure and 
magnetic pressure, and the value is related to the strength of 
the outflow. For cases considered here (s — 0.5), this critical 
value is 1. 

From the upper-right panel of Fig. (2] we see that the ac- 
cretion flow will be super-Keplerian when the toroidal mag- 
netic field is very strong (i.e., /3i > 1). This is because the 
strong magnetic stress force serves as a centripetal force. 

Figure [3] shows the case when the specific internal en- 
ergy of the outflow is different from that of the inflow. We 
take parameters a = 0.1,/ — 1,^2 = 0.01, s = 0.5 and 
^1 = 1, as before. Unlike the case shown in Fig. (2] here we 



set ^2 = 1. The dotted, solid and long-dashed lines corre- 
spond to ^3 = 0.8, 1 and 1.2, respectively. The outflow then 
plays an extra cooling role for the inflow if its specific inter- 
nal energy is higher than that of the inflow (^3 > 1), or a 
heating role otherwise. We can see from the figure that for 
the case of a weak magnetic field (/3i < 1), when ^3 increases 
from 0.8 to 1.2, the temperature (cs) and radial velocity of 
the infiow decrease while the angular velocity increases. The 
reason is as follows. From Eq. (|25|) we know ci and C3 must 
increase or decrease together. If both would increase with 
increasing ^3, from Eq. (|24|) we know that C2 must decrease, 
which is in conflict with Eq. (|27|l . The results of a strong 
magnetic field (/?i > 1) can be understood in a similar way. 

Unlike the case without outfiows (Fig. [TJ , we find from 
Figs. [2] & [3] that, when the toroidal component of the large- 
scale magnetic field is strong (e.g., /3i ~ 10), the accretion 
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Figure 3. The dynamics of tlie accretion flow wlien tlie outflow lias different specific internal energy from that of the inflow. Parameters 
are a = 0.1,/ = l,s = 0.5,(32 = 0.01 and = ^2 = 1- The dotted, solid, and long-dashed lines correspond to ^3 = 0.8,1.0, and 1.2, 
respectively. The bottom-right panel shows the ratio (C4; Eq. ([28J) of the angular momentum transport due to MRI, to that due to the 
ordered magnetic field. 



flow is super-Keplerian. This is because, for the cases of 
strong outflows (s = 0.5; M oc r) considered in Figs. [2] & 
[S] the "centrifugal force" (gradient of the gas plus magnetic 
pressure) is much weaker than in the case without outflows 
(Fig- HI); while the "centripetal force" (magnetic stress force) 
does not vary too much. 



Finally, we note that in the limits of /, /32 ^ 0, s ^ 
—0.5 and ^1,^2,^3 —> 1, i.e., the accretion flow is radiatively 
efflcient, and is devoid of outflows and magnetic flelds, our 
equations will lead to the solution with — » 0, Cs ^ 
and vk, which corresponds to a standard thin disc, as 

expected. 



5 SUMMARY 

Numerical MHD simulations show that both large-scale or- 
dered magnetic flelds and outflows exist in the inner regions 
of an ADAF. In this paper, we have investigated their influ- 
ences on the dynamics of the accretion flow in a self-similar 
approach. We assume that in addition to the stochastic com- 
ponent, the magnetic field has a strong large-scale ordered 
component. The magnetic pressure and stress forces associ- 
ated with the large-scale field thus will affect the dynamics of 
the accretion flow. We adopt the conventional a description 
to mimic the angular momentum and heating effects asso- 
ciated with the tangled magnetic field. We find that when 
the large-scale magnetic field gets stronger, the radial and 
angular velocities of the accretion fiow increase, while the 
temperature decreases. The large-scale field will be more im- 
portant than the turbulent-viscosity in transferring angular 
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momentum even when the field is moderately strong (Fig. 

We have obtained the self-similar solutions with both 
large-scale magnetic fields and outfiows. The outfiow does 
not only change the radial profile of the mass accretion rate, 
but is also able to take away angular momentum and energy 
from the inflow. We parameterize the physical properties 
of the outflow, namely its radial velocity, specific angular 
momentum and specific internal energy. We find that when 
the specific angular momentum of the outflow increases, the 
temperature of the inflow decreases, while the radial and an- 
gular velocities increase (decrease) when the magnetic field 
is weak (strong) (Fig. [2]). When the specific internal en- 
ergy of the outflow increases, the temperature and radial 
velocity of the inflow decrease, while the angular velocity 
increases (decreases) when the large-scale magnetic field is 
weak (strong) (Fig. [3]). When both outflows and large-scale 
magnetic flelds are present, we flnd that the inflow could be 
super-Keplerian if the fleld is very strong. 

The decrease in temperature of the inflow in presence of 
large-scale magnetic flelds and/or outflows could explain a 
puzzle that the predicted temperature of hot accretion flows 
is higher than that obtained from fitting the observational 
data. 
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